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Summary. The benzodiazepine receptor (peripheral) (BZRP) mainly local- 
ized on glial cells plays a role in neurosteroid synthesis, and increases with 
glial proliferation. We have recently reported a significant decrease in the 
density of BZRP labeled by [ 3 H] PK 11195 in the postmortem brain of chronic 
schizophrenics, suggesting that dysfunctions of the BZRP are involved in the 
pathophysiology of schizophrenia. We screened 11 patients with schizophre- 
nia and 10 controls, which were used in a previous postmortem study, for their 
genomic sequences of the BZRP gene in order to find DNA sequence varia- 
tions. One novel missense polymorphism (Hisl62Arg) and another previously 
reported missense mutatio n (Alal47Thr) wer e detected. An association study 
of the identified variationlTwas then performed in an extended Japanese 
sample of 304 schizophrenic patients and 369 controls. While there was an 
increased tendency in the frequency of lhe_J 6? A rp allele of schizophren ics 
compared to t hat of the controls fp = 0.0603). no statistically significant^ 
association with schizophrenia was observed in th e Alal47Thr allele (p = 
0.1016). These results do not suggest that the two missense polymorphisms 
play a major role in the genetic predisposition of schizophrenia in the 
Japanese sample. 


Keywords: Benzodiazepine receptor (peripheral), schizophrenia, single 
nucleotide polymorphism, association, missense, restriction fragment length 
polymorphism. 


♦Present address: Laboratory for Molecular Psychiatry, Brain Science Institute, 
Riken, 2-1 Hirosawa, Wako-city, Saitama 351-0198, Japan 
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Introduction 


The benzodiazepine receptor (peripheral) (BZRP) is distributed in the pe- 
ripheral tissues, such as the adrenal gland and testis, as well as in the central 
nervous system (De Souza et al., 1985; Kurumaji and Toru, 1996). The recep- 
tor plays an important role in the rate-limiting step of steroidgenesis, which 
localizes in the mitochondrial outer membrane of steroid-producing cells 
(Kreuger and Papadopoulos, 1990). The BZRP mainly localized on astroglial 
cells (Itzha et al., 1993) is also involved in steroid genesis in the brain mito- 
chondrial preparation (McCauley et al., 1995; Papadopoulos et al., 1992). The 
neurosteriods putatively elicit potent short-term allosteric modulatory effects 
on the action of GABA at the GABA a receptor (Majewska, 1992), which may 
produce an anxiolytic effect (Auta et al., 1993), and long-term genomic- 
mediated effects acting at the intracellular steroid receptor (McEwen, 1991). 

BZRP has been shown to be increased in the brain following ischaemic 
(Stephenson et al., 1995) and excitotoxic neuronal damages (Benavides et al., 
1987), and some neurodegenerative diseases such as Huntington’s disease 
(MeGmer and Reynolds, 1998; Schoemaker et al., 1983). We have reported 
that the density of BZRP labeled by [ 3 H] PK 11195 in the postmortem brains 
of chronic schizophrenic patients was decreased in three brain areas (superior 
parietal cortex, primary visual area and putamen) within 26 examined brain 
regions, and that the maximum number of binding sites (Bmax) as well as the 
dissociation equilibrium constant (Kd) were reduced in the brain areas of 
schizophrenics (Kurumaji et al., 1997). We have speculated that the reduced 
density of BZRP in the postmortem brains of schizophrenics may be involved 
in the pathophysiology of schizophrenia, associated with dysfunctions of 
neurosteroids, and that the decreases in both Bmax and Kd might be attrib- 
uted to unknown alterations in the conformation and/or structure of BZRP in 
schizophrenia (Kurumaji et al., 1997). 

The BZRP gene, localized in the 22ql3.3 band of the human genome 
(Riond. et al., 1991), consists of four exons, with exon 1 encoding only a short 
5' untranslated segment (Lin et al., 1993). Several studies have reported some 
evidence for linkage between schizophrenia and genetic markers located on 
chromosome 22q (Coon et al., 1994; Schizophrenia Collaborative Linkage 
Group for Chromosome 22, 1998). In the present study, we screened the 
coding sequence of the BZRP gene from exon 2 to exon 4 in 11 patients 
with schizophrenia and 10 controls used in the postmortem brain study 
(Kurumaji et al., 1997) to find variations resulting in a significant change in 
the BZRP functions. The two missense polymorphisms of the nucleotide 
sequence in exon 4 was detected, one of which is a novel one. Furthermore, we 
conducted an association study between schizophrenia and these missense 
polymorphisms. 


Materials and methods 

Subjects 

The schizophrenics examined in this study were 304 unrelated Japanese patients (175 
males and 129 females) with a mean ± S.D. age of 44.5 ± 12.5, including 11 patients used 
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in the previous postmortem brain study (Kurumaji et al., 1997). All patients met the 
criteria for schizophrenia of the American Psychiatric Association’s Diagnosis and Statis- 
tical Manual of Mental Disorders (DSM-III-R) (1987). The patients were receiving 
treatment at eight hospitals within 200 km of Tokyo. The control subjects were 369 
unrelated Japanese (199 males and 170 females) including 10 controls used in the previous 
study (Kurumaji et al., 1997) with a mean ± S.D. age of 45.5 ± 12.5. Among the controls, 
203 were paramedical staff members documented to be free of psychosis. The reminder 
were corporate employees who had requested annual physical examinations but had not 
been evaluated for psychiatric disorders by a psychiatrist. 

This study was approved by the Ethics Committee of the Tokyo Medical and Dental 
University. Informed consent was obtained prior to the study. 

DNA analysis 

Genomic DNAs from all subjects were prepared from peripheral whole blood cells and 
postmortem brains by the phenol extraction method or by using a DNA Extraction Kit 
(STRATAGENE). The fragments including exons 2, 3 and 4 were amplified by the touch 
down polymerase chain reaction (PCR) using the primers and PCR conditions shown 
in Table 1. The PCR reaction mixture contained a total volume of 60 pi consisting of 
lOOng genomic DNA, 200pM of each dNTP, 5pM primers, 1.5 units of TakaRa Ex Taq 
(TAKARA SHUZO) in an equipped buffer. The “touchdown” PCR was done as follows: 
After denaturing at 95°C for lmin, 15 sec at 94 °C, 30 sec at T* plus 11 minus N** °C (T* 
is listed in Table 1 and N** is cycle number) and lmin at 72°C for the first 10 cycles, and 
10 sec at 94°C, 30 sec at T* °C and lmin at 72°C in the subsequent 20 cycles. Final 
elongation was 5 min at 72°C. Direct sequencing was carried out using the Big Dye 
terminator (PE Applied Biosystems, Norwalk, CT) on an ABI 377 auto sequencer (PE 
Applied Biosystems, Norwalk, CT). 

Genotyping of the DNA sequence variants was performed by means of PCR-based 
restriction fragment length polymorphism (RFLP) assays (Table 2). A 30-pl aliquot of the 
PCR product was incubated with 2U of Nru I or 2U of Nla III according to the 
manufacturer’s recommendations (New England Biolabs, Schwalbach, Germany). Frag- 
ments were separated in a 3% NusieverSeakem agarose gel, and the bands were visualized 
by ethidium bromide staining and ultraviolet transillumination. 

Statistical procedures 

Deviation of the genotype counts from Hardy-Weinberg equilibrium were tested with 
HWE using Linkage Utility Programs (Lathrop et al., 1984). Differences between the 
patient with schizophrenia and the controls in allele frequencies and in genotype distribu- 
tions were analyzed by the x 2 test and by the Monte Carlo method using the CLUMP 
program (Sham and Curtis, 1995), respectively. 

Results 

The BZRP gene consists of four exons. Exon 1 encodes a short 5' untranslated 
segment, and exons 2, 3 and 4 code 169 amino acids of the BZRP protein. We 
screened 11 schizophrenics patients and 10 controls used in the previous 
postmortem brain study (Kurumaji et al., 1997) for sequence variations within 
the coding region (exons 2-4) using the direct sequence method. One novel 
missense variant in exon 4, which is derived from the nucleotide transition in 
codon 162 (CGT to CAT) resulting in the histidine to arginine fHisl62Ar pl 
change (Fig. 1), was detected i n threejof 11 schizophrenics a nd 4 of 10 con- 
trols. Another missense variant in exon 4, which was previously reported 
(Mirel DM, accession ID of Genome Data Bank is 9861863), was also found 
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The primers were designated from the sequences deposited in the Genbank data base (Accession nos. L21952, L21953, L21954) 
* Final annealing temperature in the touch down method (see text) 
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Table 2. Detection of DNA sequence variants in exon 4 of human benzodiazepine receptor (peripheral) 

gene 


Variant 

Primer pair 

PCR product (bp) 

Restriction enzyme 

Allele 

Fragment size (bp) 

CG/AG at 

Exon 4F & 

329 

Nla III 

CGG 

106 + 223 

codon 162 

Exon 4R-1 



CAG 

9 + 97 + 223 

G/ACG at 

Exon 4F & 

299 

Nru I 

GCG 

115 + 184 

codon 147 

Exon 4R-2 



ACG 

299 


i 


rp nnr' r^T ; 
i. \JrOV-v 


TGQCNTGGG 



M 1 2 3 


Fig. 1. Sequence and PCR-RFLP pattern of the Hisl62Arg polymorphism. The PCR 
product amplified using the Exon 4F and Exon 4R-1 primers (Table 1) was digested with 
two units of restriction enzyme Nla III. The fragments were separated on 3% agarose gels. 
Lane 1, individual homozygous for CGG (Hisl62His); lane 2, individual heterozygous for 
C G/A G (Hisl62Arg); lane 3, homozygous for CAG (Argl62Arg); lane M, 100 bp DNA 

Ladder (GIBCO BRL) 


in one of 10 control subjects, which is a change from GCG to ACG at codon 
147 with a substitution of alanine with threonine (Alal47Thr, Fig. 2). Al- 
though we further screened exon 2, exon 3 and exon 4 in an additional 40 
schizophrenic patients by sequencing, no variation in the DNA sequences was 
observed. 
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TTCCCG/CC TTC HCGPC 



300 bp 
200 bp 


100 bp 



M 1 2 


Fig. 2. Sequence and PCR-RFLP pattern of the Alal47Thr polymorphism. The PCR 
product amplified using the Exon 4F and Exon 4R-2 primers (Table 1) was digested with 
two units of restriction enzyme Nru I. The fragments were separated on 3% agarose gels. 
Lane 1, individual heterozygous for G/A CG (Alal47Thr); lane 2, homozygous for GCG 

(Alal47Ala); lane M, 100 bp DNA Ladder (GIBCO BRL) 


The genotype and allele distributions of Hisl62Arg in the schizophrenics 
and controls are shown in Table 3. There was no statistically significant 
difference between the schizophrenics and controls either in genotype (x 2 = 
4.323, d.f. = 2, p = 0.120) or allele frequencies (x 2 = 3.349, d.f. = 1, p = 0.060) 
of the Hisl62Arg variant. The frequencies of the 162Arg allele were 25.7% in 
the schizophrenic patients and 21.2% in the controls. Moreover, no statisti- 
cally significant difference between the schizophrenics and controls was ob- 
served in the distribution of the genotype (x 2 = 3.185, d.f. = 2, p = 0.103) or 
allele (x 2 = 2.543, d.f. = 1, p = 0.102) of Alal47Thr (Table 4). The frequency 
of individuals carrying 147Thr was 3.8% for the schizophrenic group and 2.2% 
for the control group. 

The distribution of genotypes in each variant did not significantly differ 
from the value expected according to the Hardy-Weinberg equilibrium: 
Hisl62Arg [controls (x 2 = 0.704, p = 0.4013), schizophrenics (x 2 = 3.269, p = 
0.0706)], Alal47Thr (controls (x 2 = 0.181, p = 0.6703), schizophrenics (x 2 = 
0.470, p = 0.4930)]. 
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Table 3. Allele frequency and genotype count for the Hisl62Arg polymorphism of the 
benzodiazepine receptor (peripheral) gene in the controls and schizophrenics 


Phenotype 

Genotype count 



P 

n 

162 Arg/ Arg 

162His/Arg 

162His/His 

Controls 

369 

14 

129 

226 

0.120* 

Schizophrenics 

304 

14 

128 

162 


■ 

Allele count 





n 

162 Arg 


162His 


Controls 

738 

157 


581 

0.0603 b 

Schizophrenics 


156 


452 



“X 2 = 4.323, d.f. = 2, b x 2 = 3.349, d.f. = 1 


Table 4. Allele frequency and genotype count for the Alal47Thr polymorphism of the 
benzodiazepine receptor (peripheral) gene in the controls and schizophrenics 


Phenotype 

Genotype count 



P 

n 

147Thr/Thr 

147Ala/Thr 

147 Ala/ Ala 

Controls 

369 

0 

16 . 

353 

0.103 a 

Schizophrenics 

304 

0 

23 

281 



Allele count 





n 

147Thr 


147Ala 


Controls 

738 

16 


722 

0.1016 b 

Schizophrenics 

608 

23 


585 



a X 2 = 3.185, d.f. = 2, b x 2 = 2.543, d.f. = 1 


Discussion 

The present study revealed a novel missense polymorphism (Hisl62Arg), and 
confirmed the previously reported missense variation (Alal47Thr) in exon 4 of 
the BZRP gene. Although there was a tendency of increased frequency of the 
162 Arg allele in the schizophrenic patients compared to that of the controls, no 
statistically significant difference between the schizophrenics and controls was 
detected in the distribution of the genotype or allele of both variation sites. 
These results suggest that the polymorphisms do not play a major role in the 
genetic predisposition of schizophrenia in the Japanese cohort. 
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The BZRP shows an ubiquitous distribution in the central nervous system 
and in the peripheral tissues (De Souza et al., 1985; Kurumaji et al., 1996). The 
receptor localizes in the mitochondrial outer membrane of the steroid- 
producing cells (Anholt et al., 1986; Antkiewicz-Michaluk et al., 1988), and 
plays an important role in the translocation of cholesterol from the outer 
to inner mitochondrial membranes, a rate-limiting step in steroidgenesis 
(Kreuger and Papadopoulos, 1990). The cDNA of the BZRP gene encodes a 
18kDa protein which has five potential transmembrane regions (Sprengel et 
al., 1989). It is suggested that the 147 Ala situates in the fifth putative trans- 
membrane region, and that the 162His localizes in the carboxyl-terminal 
sequence. The 147 Ala is conserved in BZRP of human, bovine and rodents, 
while the 162His is replaced by Arg in the receptor of bovine and rodents 
(Farges et al., 1994; Lin et al., 1993). Moreover, it is proposed that BZRP may 
conform a multimeric receptor complex composed of multiple 18kDa protein 
subunits and other proteins such as the voltage-dependent anion channel 
(Joseph-Liauzun et al., 1997; Papadopoulos, 1998). 

The BZRP can be labeled with isoquinoline carboxamides such as PK 
11195 and with benzodiazepines such as Ro5-4864, although the binding 
domains of the compounds are overlapping but not identical sites on the 
BZRP (Farges et al., 1993; Rao and Butterworth, 1997). The deletion of the 
residue 15-35 sequence of the amino-terminal of the BZRP resulted in 
the complete loss of- binding of both [ 3 H] PK 11195 and [ 3 H] Ro5-4864, 
whereas the deletion of the last 14 amino acids of the receptor at its carboxyl 
terminus abolished the binding of Ro5-4864, with no effect on the binding of 
PK 11195. A particular residue of the sequence, 154Val, was crucial for the 
binding of the benzodiazepine ligand, but not of the isoquinoline ligand 
(Farges et al., 1994). In addition, the Hisl62Arg and Alal47Thr variations 
determined in the subjects used in the postmortem brain study were not 
correlated to the changes in the binding parameters of [ 3 H] PK 11195 in the 
brains (Kurumaji et al., 1997) (data not shown). 

In conclusion, there were no associations between the two missense poly- 
morphisms of the BZRP gene and schizophrenia in the Japanese cohort. 
Taking into account of the heterogeneity of schizophrenia, it may be impor- 
tant to investigate associations between these polymorphisms and schizophre- 
nia in other populations. Moreover, it remains to be clarified whether or not 
the variations in the amino acid sequence of BZRP result in an altered 
interaction between BZRP and other mitochondrial membrane proteins and 
in a change of the physiological functions of the receptor complex. A search 
for other variations in the BZRP gene is also needed, particularly in the 
regulatory region, which can produce a change in the density of BZRP in the 
schizophrenic brains. 
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Tam, 

You had requested the public availability date for a sequence with accession # 9861863. I was unable to find 
such a record. 

I searched the protein, genome, and nucleotide databases at NCBI (via Entrez-PubMed). There were no hits in 
the protein & genome databases. The nucleotide database did have a match to that accession number, but it does 
not appear to be the correct sequence. The author "Mirel DM" does not appear in the bibliographic information, 
and the sequence is from wheat. I've attached the record below so you can review it for yourself. 
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The document below contains screen shots of other search attempts. 

Page 1 is the result of searching '9861863' at the DDBJ website (no matches) 

Page 2 is the result of searching '9861863' at the Genomic Data Bank (no matches) 

Page 3 is the result of searching '9861863' AND 'Mirel' in Google. There was only one hit, which appears to be 
the article you already have. 
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I believe the accession number citation is not correct. 

Please let me know if you have any questions. 
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dbEST Id: 

EST name: 
GenBank Acc : 
GenBank gi : 

CLONE INFO 
Clone Id: 

DNA type : 

PRIMERS 
Sequencing : 
PolyA Tail : 

SEQUENCE 


Entry Created: 
Last Updated: 

COMMENTS 


LIBRARY 
dbEST lib id: 
Lib Name: 
Organism: 
Cultivar : 

Tissue type: 
Develop . stage : 
Lab host : 
Vector: 

R. Site 1: 

R . Site 2 : 
Description: 


SUBMITTER 

Name: 

Institution : 

Address : 

Tel : 

Fax : 


5813432 

WHE14 05 - 14 08_J03_J03ZS 

BE604593 

9861863 


WHE1405-1408_J03_J03 

CDNA 


Stratagene SK primer 
Unknown 


TGTGTGTACCACGACCACCCCTTTTATCTCATCCCATCCAAGCTAGACCAGTCCCGCACT 
GTGCTGTGGTGCTCCATCTAGCCAGCTTTTTTGCTCAGCTTAGCTTCCAAATCTGCTCCA 
CAGCTCCGGCGGCACAAGAAGAAAATATCCGTGCGTCCGCGAGAGGAAGGGAGATCGTCC 
ATGGAGAAGGACAACCACAAGGAGGGGGAGGGGG CGACGG CGTGGATGACGGTG CCGG CG 
TTCGGGGACTGGGACATGAAGAACGGCGCCATGCCGGACTACTCCATGGACTTCTCCAAG 
ATCCGGGAGATG CG CAAG CAGAACAAGAAGGAGCTC 

Aug 21 2000 
Aug 21 2000 


Sequence have been trimmed to remove vector sequence and low 
quality sequence with phred score less than 20 

5596 

Wheat drought stressed leaf cDNA library 
Triticum aestivum 
TAM W-101 
Leaf 

Full tillering stage 
E. coli SOLR 

Lambda Uni -ZAP XR, excised phagemid 

EcoRI 

Xhol 

Plants were given a gradual stress down to 65% and 78% RWC 
at Texas Tech University (D. Zhang in HT Nguyen lab) . Total 
RNA and poly (A) RNA were prepared, a cDNA library was made, 
and the cDNA clones were in vivo excised to give pBluescript 
phagemids in the TJ Close lab (Choi, Close), at the 
University of California, Riverside. Plasmid DNA 
preparations and DNA sequencing were performed in the OD 
Anderson lab (all other authors) . 
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01 in Anderson 

US Department of Agriculture, Agriculture Research Service, 

Pacific West Area, Western Regional Research Center 

800 Buchanan Street, Albany, CA 94710, USA 

5105595773 

5105595818 



E-mail : 

oandersn@pw . usda . gov 
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Title: 

The structure and function of the expressed portion of the 
wheat genomes - Drought stressed leaf cDNA library 

Authors : 

Anderson, 0. D. , Chao,S., Choi,D.W., Close, T.J., Han,P.S., 
Hsia,C.C., Kang,Y., Lazo,G.R., Miller, R., Nguyen H.T., 
Rausch, C.J., Seaton, C.L., Tong,J.C., Zhang, D. 

Year: 

2000 

Status : 

Unpublished 
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